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ABSTRACT. We investigate some qualitative behavior of the solutions of the difference e- 
k k 

quation £n+1 = a+ DS) a(2441)2n—(2i41) /b¥n—e+ DS B(2i)n—(2i) Where the the initial condi- 
i=0 i=0 


tions £-r,%-r+1,--., Zo are arbitrary positive real numbers such that r = max {¢,k} where 


i,r € {0,1,...} and a, Q(2i+1); Bai) are positive constants. 


1. INTRODUCTION 


In this paper we deal with some properties of the solutions of the difference equation 
k 
a Q(2i+1)Tn—(2i+1) 


(1.1) Ae e 





7 s w20 Ti 2s 
bine + D7 B(2i)En—(28) 
i=0 


where the the initial conditions x_,,2%_,41,...,% are arbitrary positive real numbers such 
that r = max {¢,k} where i,r € {0,1,...} and a, ai41), B27) are positive constants. There is 
a class of nonlinear difference equations, known as the rational difference equations, each of 
which consists of the ratio of two polynomials in the sequence terms in the same form. there 
has been a lot of work concenring the global asympototic of solutions of rational difference 
equations [2], [3], [4], [5], [6], [7], [8], [10] and [11]. 

Many reseaches have investigated the behavior of the solution of difference equation for 
example: 

Amleh et al. [1] has studied the global stability, boundedness and the periodic character 


of solutions of the equation 
In-1 





Ins =Q + ig 
n 


Our aim in this paper is to extend and generalize the work in [1]. That is, we will investigate 
the global behavior of (1.1) including the asymptotical stability of equilibrium points, the 
existence of period two solution and investigate the oscillation property of the recursive 
sequence of Eq. (1.1). 

Now we recall some well-known results, which will be useful in the investigation of (1.1) 
and which are given in [9]. 
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Let I be an interval of real numbers and let 
F: ee 
where F is a continuous function. Consider the difference equation 
(1.2) tniii E (an Taams nka 0 Plea 
with the initial condition £k, £-k41; -£o € I. 


Definition 1. (Equilibrium Point) 


A point 7 € J is called an equilibrium point of Eq. (1.2) if 





T= f(T Test): 
That is, £n = 7 for n > 0, is a solution of Eq. (1.2), or equivalently, © is a fixsed point of f. 


Definition 2. (Stability) 


Let % € (0,00) be an equilibrium point of Eq. (1.2). Then 
(i) An equilibrium point 7 of Eq. (1.2) is called stable if for every £ > 0 there exists ô > 0 
such that, if £r, £-r41, ---, £o € (0,00) with |x_, — T|+|£-r41 — F]+...+|v9 — Z| < 4, 
then |x, — T| < € for all n > —r. 
(ii) An equilibrium point z of Eq. (1.2) is called locally asymptotically stable if © is 
locally stable and there exists y > 0 such that, if v_,,@_,41,...,%0 € (0,00) with 
[£r — T| + |£-r41 — T| +... + [£0 — T| < y, then 


lintr = @. 


(iii) An equilibrium point 7 of Eq. (1.2) is called a global attractor if for every £-r, G41, ...,Z0 € 
(0, 00) we have 


lim £n =T. 
noo 


(iv) An equilibrium point % of Eq. (1.2) is called globally asymptotically stable if z is 
locally stable and a global attractor. 
(v) An equilibrium point % of Eq. (1.2) is called unstable if 7 is not locally stable. 


Definition 3. (Permanence) 


Eq. (1.2) is called permanent if there exists numbers m and M with 0 < m < M < œ 
such that for any initial conditions £-p, %-r+1, ---, £o E (0,00) there exists a positive integer 
N which depends on the initial conditions such that 


m< zn ISM foralln > —N. 


Definition 4. (Oscillation) 
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A sequence {£n} _, is called nonoscillatory about the point 7 if there is exists N > —k 
such that either x, > for all n > N or x, < F for all n > N. Otherwise {x,,}°°_, is called 
oscillatory about 7T. 


Definition 5. (Periodicity) 


oo 
n=—-r 


is said to be periodic with prime period p if p is the smallest positive 


A sequence {z,,} is said to be periodic with period p if ,4, = x, for all n > —r. 


OO 
n=-r 


A sequence {£n} 
integer having this property. 

The linearized equation of Eq. (1.2) about the equilibrium point % is defined by the 
equation 


k 
(1.3) Zn+1 = S Dace 
i=0 


where 
OF(T,T,..., T, ) 
OLn—i 





Pi = ; i = 0,1,..., k. 


The characteristic equation associated with Eq. (1.3) is 
(1.4) AFHL L poA? — p AFT! =n — ppi À — pk = O. 


Theorem 1.1. [9]. Assume that F is a C! — function and let T be an equilibrium point of 
Eq. (1.2). Then the following statements are true: 
(i) If all roots of Eq. (1.4) lie in the open unit disk |A| < 1, then he equilibrium point z 
is locally asymptotically stable. 
(ii) If at least one root of Eq. (1.4) has absolute value greater than one, then the equilib- 
rium point © is unstable. 
(iii) If all roots of Eq. (1.4) have absolute value greater than one, then the equilibrium 
point T is a source. 


Theorem 1.2. [12] Assume that p; € R,i = 1,2, ..., k. Then 


k 
X` Ir: <1, 
i=l 


is a sufficient condition for the asymptoticcally stable of Eq. (1.5) 


(1.5) Yntk + PiYn+k-1 t+ + PkYn = 0, n = 0,1,.... 
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2. LOCAL STABILITY OF THE EQUILIBRIUM POINT OF EQ.(1.1) 


In this section we investigate the local stability character of the solutions of Eq. 


Eq. (1.1) has a unique nonzero equilibrium point 











k 
D aeiy) 
ah i=0 
T=ac 7 , 
bx + yy Bait 
i=0 
then 
k 
De Q(2i+1) 
T=0a+ = 
b+ 3 Bi 
Let 
k k 
G=} aci and Q= >) Bay. 
i=0 i=0 
Then, we get 
T=at z 
b+Q 
Let f : (0,00)*t' — (0,00) be a function defined by 
k 
2 O(2i+1)U (9441) 
(2.1) f (v, uo, Uy, Ue) = a+ = 


E à 
bu + 5 B eguei) 
i=0 


Therefore it follows that 














k 
OF (V, uo, U1, Uk) E 2, Cit) 
Oui = k j 
(2i+1) [be +5 Benue] 
i=0 
k k 
Of (V, Uo, U1, ---5 Uk) _ T> Pe e aitosa) 
ua) k 2 , 
e +208 oat 
i=0 


and 


k 
—b i , 
Of (v, U0, U1, «5 Uk) = (sae 1a) 


Ov a k z 
i=0 
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Then we see that 











Of (TT; , T) = G Poi 
ugi)  a(b+Q)+G 2+1; 
Of (T,F, 2,7) _ -GQ ae 
Dua BrO 
and 
OF Tasa bG =P 
dvs (+ Q)(a(B4+Q)+G) M 


Then the linearized equation of (1.1) about Z is 


k 
(2.2) Zn+1 = N Pini: 
i=0 


Theorem 2.1. Assume that 
G<a(b+Q). 


Than the equilibrium point of Eq. (1.1) is locally stable. 


Proof. It is follows by Theorem(1.2) that, Eq. (2.2) is locally stable if 





pel +. + [pg] + -+ [pil + [pol < 1. 


That is 


e CERO TCCES Exo] hal (CECT 
a(b+Q)+G] |(6+ Q)(a(b+Q)+G)| |(b+Q) (a(b+ Q) +G) 





<1, 





this implies that 
bG + G (b+ Q) + GQ < (b+ Q) (a(b+ Q) +G). 
Thus 
G<a(b+Q). 











Hence, the proof is completed. 





Example 2.1. Consider the difference equation 


0.624 + 0.2£n-3 


n = 0.54 ; 
Fn E E N 








where k = 1,4 = 2,a = 0.5,b = 2,a1 = 0.6, @3 = 0.2, 8o = 9, 8. = 0.5. Figure(2.1), shows 
that the equilibrium point of Eq. (1.1) has locally stable, with initial data £n-3 = 2.2, £tn-2 = 





6 ELABBASY, BARSOUM AND EL-BIATY 


0.2, 2n-1 = 3.2, £o = 0.6. 


x(i+1)=0.5+(0.6°x(i-1)+0.2*x(i-3)/(2*x(i-2)+9*x(i)+0.5*x(i-2))) 























Figure 2.1 


3. GLOBAL STABILITY OF EQ. (1.1) 


Our aim in this section we investigate the global asymptotic stability of Eq. (1.1). 

Considar of the Eq. (1.2) we give the following theorem which is a minor modification of 
Theorem (A.0.6) in [9] and which will be useful for the investigation of the global attractivity 
of solutions of Eq. (1.1). 


Theorem 3.1. Let [a,b] be an interval of real numbers and assume that 
Fed = [a, 8] 


is a continuous function satisfying the following properties: 
(a) f (@1,%,..., p41) is non-decreasing in the ugi+41 for each v and uy; in [a,b] and non- 
increasing in the v and ug; for each uz +1 € |a, b] ; 
(b) If (m, M) € [a,b] x [a,b] is a solution of the system 
m= f(M,m,M,m, M, m, ..., M,m, M), 
and 
M = f (m, M,m, M,m, M,....m,M,m), 
then 
m= M. 
Then Eq. (1.2) has a unique equilibrium T € [a,b] and every solution of Eq. (1.2) 


converges to T. 


Proof. Set 


mo=a and Mj=b 
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and for each i = 0,1, 2, ...Set 


m; = f (M; 1s Mi—1, Mi—1, Mi—1, Mi—1, Mi b Mii, Mi-1, Mi-1), 





and 
M; = f (m; 1; Miya, Maas Mi-1, Mi 1) e Mi—1, Mi—1, Mi-1) - 





Now observe that for each i > 0, 


a = mı < Mmg < Ms <... < mi <... < Mi < ... < Ma < Mo < Mo = b, 





and 
M; Sı SM, fo q>(k+1)i+1. 
Set 
m= lim m; and M = lim M;. 
ioo i=oo 
Then 


M > lim sup 2; > lim inf z; > m, 


1—00 1—00 


and by the continuity of f, 

m= f(M,m,M,m, M, m, ..., M,m, M), 
and 

M = f (m, M,m, M,m, M, ...,m, M, m). 


In view of (b), 
m= M, 











from which the result follows. 





Theorem 3.2. Ifa(b+ Q) 4G, then the equilibrium point T of Eq. (1.1) is global attractor. 


Proof. Let f : (0,00)**! — (0,00) be a function defined by 


k 
D Q(2i+1)U(2i+1) 





, 


k 
bu + > B (iy 422) 
i=0 


then we can see that the function f (v, uo, u1, ..., Uz) is decreasing in the v, ug; and increasing 
IN Uzi+1- 
Suppose that (m, M) is a solution of the system 


je = f(M,m,M,m, My ..., M,m, M 


and 
M = f (m, M,m, M,m, M,...,m,M,m). 
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Then from Eq. (2.1), we see that 











k k 
5 Q (2i+1)M > Q (i+1) M 
i i=0 n i= 
m = aq k , M=a4 k ? 
bM + >) BaiyM bm +d) bem 
i=0 i=0 
at Gm M T GM 
m = a+ ——— =a+ ———, 
bM+QM’ bm + Qm 


bMm+QMm abM +aQM+Gm, 
bMm+QMm = abm+aQm+GM, 





than 
ab(M —m)+aQ(M —m) -—G(M —m) =0. 
Thus 
m= M. 


It follows by Theorem(3.1) that % is a global attractor of Eq. (1.1) and then the proof is 
complete. 














4. OSCILLATION OF EQ. (1.1) 


Theorem 4.1. If k — odd, — even, then Eq. (1.1) has an oscillatory solution 


Proof. firstly if k > € we assume that, 





(4.1) Dp EA bE ee LL, and U_k41,U—-k43, +0, LQ < 2. 
So, 


k 
D Q (21-41) —(2i+1) 
a = a + ——— 5 
bx_¢ + 5 BT- 
i=0 





then 





and 
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So, we have 


k 
Ly QO (2141) U—(2i41)41 
i=0 





Ly = at ’ 


k 
bze +) beyt-ei+ 
= 


SO, 





then, 


Secondily assume that, 


Lok opi te <T, and eps k43; To > T. 

















So, 
k 
>D Q (2+1) T- (2i+1) 
zı = a + —— ; 
bx_¢ + 2 aie 
then, 
k 
ss Q(2i41)U 
t<ac = k ; 
bT + 2 Bent 
and 
k 
X Q(2i+1) 
t,<ac =e i = 
b+ 2 Bi 
So, we have 
k 
oD O (2141) —(2i41)41 
z2 = a + — : 


k 
bze +) bent- 
i=0 
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SO, 
k 
> Q (2i+1)T 
Lg > a + ~ 7 f 
bT + > begt 
i=0 
then, 
k 
> Q(2i+1) 
£o >a + = 7 = 
b+ > Pin 


One camproceed in prove manwer to show that £3 < % and z4 > ¥ and soon. Hence, the 











proof is completed. 





Example 4.1. Consider the difference equation 

0.0017,_1 + 0.52,_3 
0.1%,-2 + 0.82, + 0.0012,’ 
where k = 1, £ = 2,a > 05,6 = U1, o¢ = 0.001, a3: = 0.5,,8, = 0.8785. = 0.001. figure(4.1) 
shows that the solution of Eq. (1.1) oscillates about T = 1.0559. Where the initial data 
satisfies condition(4.1) of Theorem(4.1) x_3 = 0.4, x2 = 2.1, £—1 = 1.3, zo = 0.2. (see Table 
41) 





Ln+1 = 0.54 
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1 | 0.4000 ||| 21 | 1.2524 ||| 41 | 1.2785 ||| 61 | 1.4099 81 | 1.6305 
2 | 2.1000 ||| 22 | 0.9259 ||| 42 | 0.8926 ||| 62 | 0.8348 82 | 0.7671 
3 | 1.3000 ||| 23 | 1.1709 ||| 43 | 1.2902 ||| 69 | 1.4271 83 | 1.6598 
4 | 0.2000 ||| 24 | 0.9486 ||| 44 | 0.8880 ||| 64 | 0.8283 84 | 0.7601 
5 | 1.0388 ||| 25 | 1.2358 ||| 45 | 1.2998 ||| 65 | 1.4451 85 | 1.6908 
6 | 1.5913 ||| 26| 0.9192 ||| 46 | 0.8823 ||| 66 | 0.8218 86 | 0.7531 
7 | 1.0037 ||| 27| 1.2056 ||| 47| 1.3122 ||| 67| 1.4642 87 | 1.7239 
8 | 0.6117 ||| 28 | 0.9363 ||| 48 | 0.8768 ||| 68 | 0.8152 88 | 0.7460 
9 | 1.2999 || | 29 | 1.2951 ||| 49 | 1.3234 69 | 1.4843 89 | 1.7581 
10| 1.1977 ||| 30 | 0.9150 ||| 50 | 0.8712 ||| 70 | 0.8085 90 | 0.7990 
11 | 0.9938 ||| 31 | 1.2805 ||| 51 | 1.3366 ||| 71 | 1.5055 Of | 1.7947 
12 | 0.8315 ||| 32 | 0.9230 ||| 52 | 0.8653 ||| 72 | 0.8018 92 | 0.7319 
18 | 1.3277 ||| 33 | 1.2446 ||| 53 | 1.3494 13 | 1.5279 93 | 1.89339 
14 | 1.0160 ||| 34 | 0.9094 ||| 54 | 0.8595 ||| 74 | 0.7950 94 | 0.7249 
15| 1.0551 ||| 35 | 1.2508 ||| 55| 1.3636 ||| 75 | 1.5515 95 | 1.8742 
16 | 0.9260 ||| 36| 0.9106 ||| 56 | 0.8534 76 | 0.7881 96 | 0.7180 
17| 1.2881 ||| 37| 1.2598 ||| 57| 1.3781 71| 1.5764 97 | 1.9173 
18 | 0.9477 ||| 38| 0.9018 ||| 58| 0.8474 78 | 0.7811 98 | 0.7111 
19| 1.1207||| 39| 1.2701 ||| 59| 1.3937 ||| 79| 1.6027 99 | 1.9630 
20 | 0.9519 ||| 40 | 0.8991 60 | 0.8411 ||| 80 | 0.7742 ||| 100 | 0.7042 
Table 4.1 


x(n) 


x(i+1)=0.5+(0.001*x(i-1)+0.5*x (i-3)/(0.1*x(i-2)+0.8*x (i)+0.001*x (i-2))) 
2.27 









































Figure 4.1 
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5. PERIODIC SOLUTIONS 
In this section we investigate the periodic character of the positive solutions of Eq. (1.1). 
Theorem 5.1. Eq. (1.1) has positive prime priod two solution only If 
(5.1) €—odd and (a—)(a8 — aa + b) > 4aap. 


Proof. Assume that there exists a prime period-two solution 























-3 P, q, P, q, sone 
of (1.1). Let £n = ¢,%n41 = P. Since l — odd, we have x,_¢ = p. Thus, from Eq. (1.1), we 
get 

EPR Qp + A3p + ... + Q(2i+1)P 
bp + Bog + Bog +...4 bena’ 
and 
PEE ang + agg +... + Q(2i+1)4 
bq + Bop + Bop +... + Bap’ 
Let 
Qi + 03+... + Q aint) O 
and 
By + Ba++ bay =À 

Than 

E YP 

me RE 

and 

_ yq 

q=at Sonor. 

Than 
(5.2) bp? + ôpq = abp + aôq + yp, 
and 
(5.3) bq? + ôpq = abg + adp+ yq. 


Subtracting (5.2) from (5.3) gives 
b (p° — q’) = (ab — aô + y) (p — q) . 
Since p Æ q, we have 
ab — aô + 
(5.4) ptq=—— L 


Also, since p and q are positive, (ab — ad + y) shuold be positive. Again, adding (5.2) and 
(5.3) yields 


(5.5) b(p? + q?) + 26pq = (ab + ad + y) (p + q). 
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It follows by (5.4), (5.5) and the relation 














p+q@=(p+q)’—2pq, VYp,qER, 


that 


ad(ab — ad + 7) 
b(5 — b) 





(5.6) pq = 


Assume that p and q are two distinct real roots of the quadratic equation 


_ab—adt+y, _ ad(ab—ad+y) _ 








2 
f b i b(5 — b) = 

and so 
ab —ad +y * _ 4ad(ab— 05+) a 
b b(d — b) 


which is equivalent to 


(6 — b)(ab — ad + y) > 4abô. 


Thus, the proof is completed. 


Example 5.1. Consider the difference equation 


ATn—1 
nt1 = 0.125 
ER $ Tn-1 = 2X5 





13 














where k = 0, €—odd, a = 0.125,b = 1,a, = 4, bo = 2. Figure(5.1), shows that Eq. (1.1) which 
is periodic with period two. Where the initial data satisfies condition(5.1) of Theorem(5.1) 
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v4 = 0.1, zo = 0.3. (seeTable 5.1) 




















































































































































































































































































































































































































1 | 0.1000 ||| 17| 0.2928 ||| 33 | 0.2690 | | | 49 | 0.2686 ||| 65 | 0.2686 
2 | 0.8000 ||| 18 | 3.5545 ||| 34 | 3.6055 ||| 50 | 3.6064 66 | 3.6064 
3 | 0.3000 ||| 19| 0.2832 ||| 35| 0.2689 ||| 51 | 0.2686 ||| 67| 0.2686 
4 | 1.4583 ||| 20 | 3.5751 ||| 26 | 38.6059 | || 52 | 3.6064 68 | 3.6064 
5 | 0.4981 ||| 21| 0.2774 ||| 37| 0.2688 ||| 53 | 0.2686 ||| 69 | 0.2686 
6 | 2.5016 ||| 22| 3.5876 ||| 38 | 3.6061 | || 54 | 3.6064 70 | 3.6064 
7 | 0.4871 ||| 23 | 0.2739 ||| 39 | 0.2687 ||| 55 | 0.2686 ||| 71 | 0.2686 
8 | 3.0038 ||| 24 | 3.5952 ||| 40 | 3.6062 ||| 56 | 3.6064 72 | 3.6064 
9 | 0.4250 ||| 25 | 0.2718 ||| 41 | 0.2687) || 57| 0.2686 ||| 73 | 0.2686 
10 | 3.2427 ||| 26 | 3.5997 || | 42 | 3.6063 ||| 58 | 3.6064 74 | 3.6064 
11 | 0.3710 ||| 27| 0.2705 ||| 43 | 0.2687) || 59 | 0.2686 ||| 75 | 0.2686 
12 | 3.3801 ||| 28 | 3.6024 ||) 44 | 3.6063 ||| 60 | 3.6064 76 | 3.6064 
13 | 0.3831 ||| 29 | 0.2697 ||| 45 | 0.2686 ||| 61 | 0.2686 ||| 77| 0.2686 
14 | 3.4664 ||| 30 | 3.6040 || | 46 | 3.6063 ||| 62 | 3.6064 78 | 3.6064 
15 | 0.3084 ||| 31 | 0.2693 ||| 47| 0.2686 ||| 63 | 0.2686 ||| 79 | 0.2686 
16| 3.5208 ||| 32 | 3.6049 ||| 48| 3.6064 |64 3.6064 ||| 80] 3.6064 
Table 5.1 

, xe ta. 125+ (44-1420) 

3.5} : 

| | 

25} 

15} 

0.5 

oll 

0 0 20 30 40 50 60 70 80 90 100 

Figure 5.1 


Remark 5.1. Note that the special cases of Eq. (1.1) have been studied in [1] when k = 
0,6 0, a, 1, Bo 1, Q2i41 = 0, Bo; = 0,2 > 1. 
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